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We present a comprehensive overview of vortex pinning in single crystals of the isovalently sub- 
stituted iron-based superconductor BaFe2(Asi_ cc P cc )2, a material that qualifies as an archetypical 
clean superconductor, containing only sparse strong point-like pins [in the sense of C.J. van der 
Beek et al., Phys. Rev. B 66, 024523 (2002)]. Widely varying critical current values for nominally 
similar compositions show that flux pinning is of extrinsic origin. Vortex configurations, imaged 
using the Bitter decoration method, show less density fluctuations than those previously observed 
in charge-doped Ba(Fei_ x Co x )2As2 single crystals. Analysis reveals that the pinning force and - 
energy distributions depend on the P-content x. However, they are always much narrower than in 
Ba(Fei_ cc Co cc )2As2, a result that is attributed to the weaker temperature dependence of the super- 
fluid density on approaching T c in BaFe2(Asi_ x P x )2. Critical current density measurements and 
pinning force distributions independently yield a mean distance between effective pinning centers 
C ~ 90 nm, increasing with increasing P-content x. This evolution can be understood as being 
the consequence of the P-dependence of the London penetration depth. Further salient features are 
a wide vortex free "Meissner belt", observed at the edge of overdoped crystals, and characteristic 
chain-like vortex arrangements, observed at all levels of P-substitution. 

PACS numbers: 74.25. Op,74.25.Sv,74.25.Wx,74.62.En,74.70.Xa 



I. INTRODUCTION 

Recent vortex imaging studies performed on iron 
pnictide superconductors show evidence for nanoscale 
inhomogeneity^ being at the origin of the low-field 
critical current density and the highly disordered 
vortex structures in these materials PEU Notably, in 
Ba(Fei_ x Co :E )2As2, the critical current density j c and 
vortex distributions imaged by Bitter decoration could 
be consistently analyzed, provided that spatial hetero- 
geneity, on a scale of several dozen nm, both of the 



critical temperature T c and the vortex line energy 
is taken to be responsible for flux pinningP At higher 
magnetic fields, of the order of several tenths of Tesla, 
nano -scale heterogeneities are inefficient in pinning flux 
lines . 1 I 13 I 14 I The critical current density is then most 
likely determined^ by the scattering of quasiparticles 
in the vortex cores associated with the presence of 
atomic-size defect s in the crystalJ^EID leading to weak 
collective pinning ft 6 * 17 * A good candid ate for t hese de- 
fects are the dopant atoms themselves! 1 * 2 * 13 * ^ The na- 
ture of the dopant atoms is essential for this media- 
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nism; charged defects lead to different scattering than 
uncharged defects.^ This weak collective pinning contri- 
bution to the critical current density manifests itself as 
a plateau-like behavior in a j c (B) plot. It is present in 
all charge-doped iron-based superconductors, as well as 
in BatFei-^Ru^AssP 

On the other hand, in isovalently substituted 
BaFe2(Asi_ ;E P :c )2 there is no indication of weak collec- 
tive pinning,^ which qualifies the material as "clean" 
with respect to charge-doped iron-based superconduc- 
tors. Given recent claims^^ 2 -! that in certain iron-based 
superconductors, the vortex configuration is more or- 
dered than what was hitherto observed SHUH it i s inter- 
esting to see whether the absence of weak collective pin- 
ning has any impact on the spatial configuration of vor- 
tices. From a magnetic force microscopy (MFM) study 
at magnetic fields up to 100 Oe, Yang et al. claim 
that vortex configurations in hole-doped single crystalline 
Bai_ ;E K ;E Fe2As2 (with x = 0.28 and x — 0.4) are more 
ordered than in, e.g. Ba(Fei_ £C Co £C )2As2P Using an anal- 
ysis method similar to that of Ref. 2, they report pinning 
forces that are one order of magnitude smaller on av- 
erage. They also claimed the observation of local tri- 
angular vortex order in the optimally doped material, 
even though closer scrutiny (see section HIE below) re- 
veals the vortex ensembles to be no more ordered than 
those in Ba(Feo.925Coo. 075)2 As2- 2 Finally, the author^ 
reported on the observation of remarkable vortex chains, 
both in underdoped (x = 0.28) and optimally doped 
(x = 0.4) Bai_ :E K iC Fe2As2. The presence of these was 
attributed to vortex pinning by the twin boundaries aris- 
ing from the orthorhombic structure, at least in the 
underdoped material. Furthermore, neutron scattering 
experiments on the vortex lattice in isovalently substi- 
tuted BaFe2(Asi_ ;E P ;E )2 were performed at T = 2 K by 
Kawano-Furukawa et al. No vortex Bragg peaks were 
found for the optimally substituted compound. However, 
after annealing the samples at 500° C, a distorted triangu- 
lar vortex lattice was observed; this became more ordered 
as the applied magnetic field was increased from 0.7 to 7 
tJHI These results suggest that the disorder responsible 
for pinning in BaFe2(Asi_ ;E P ;E )2 is extrinsic in nature. 

In this work we present and analyze sustainable 
current density measurements, magneto-optical imag- 
ing, and Bitter decoration experiments performed on 
BaFe2(Asi_ x P :E )2 single crystals with different x. In 
contrast to some reports p9 we find strong pinning, pre- 
sumably by nanoscale heterogeneity, as the only ob- 
served pinning contribution in fields up to 5 T. As in 
Ba(Fei_ x Co :E )2As2 and other materials, 6 there is no ev- 
idence for any extended triangular order in the vortex 
ensemble; thus, the strong pinning contribution in itself 
suffices to generate the extreme disorder of the vortex 
ensemble. The spatial configuration of vortices in iso- 
valently substituted BaFe2(Asi_ x P :c )2 does not present 
large vortex density fluctuations such as observed in 
charge-doped Ba(Fei_ :E Co :E )2As2 single crystals,^ a fact 
that is attributed to the different temperature depen- 



dences of the superfluid density in the two materials. 
The quantitative analysis of the vortex configurations 
in terms of the pinning energy confirms that pinning 
disorder is somewhat less effective in BaFe2(Asi_ ;E P ;E )2 
than in Ba(Fei_ x Co :E )2As2, and that it depends on the 
P-content x. Analysis of the pinning energies, pinning 
forces, and the critical current density as function of P- 
content yields consistent estimates of the effective pin 
density. This clearly decreases upon increasing the P- 
content, a behavior that tracks the composition depen- 
dence of the scattering rate in the normal state. The 
main features of our results can be understood in terms 
of the composition dependence of the vortex line energy, 
implying that local variations of the superfluid density 
is a good candidate for the origin of the vortex pinning. 
The local variation of the dopant atom density as well 
as the variation of the sustainable current density with 
composition argues against any possible spatially phase- 
separated superconducting and anti-ferromagnetic states 
of the material as being at the origin of pinning. 



II. EXPERIMENTAL DETAILS 

Experiments have been performed on 
BaFe2(Asi_ ;E P cc )2 single crystals grown by the self- flux 
method,^ and characterized using Energy Dispersive 
X-Ray Spectroscopy (EDX) and EDX mapping in a 
scanning electron microscope (SEM). Crystals with 
manifest chemical heterogeneity were discarded from 
further study. The crystals described below present 
no impurity phases, within the experimental limits of 
accuracy < 1%. 

Magnetic flux penetration in crystals with different 
substitution levels (x = 0.27 — 0.49) was charac terized 
using the magneto-optical imaging (MOI) method^III] 
before further experiments. The MOI technique notably 
allows one to discard samples with macroscopic defects, 
and also, to extract calibrated flux density profiles. The 
sustainable current density j s for magnetic fields up to 
fioH a = 500 G (50 mT) was obtained from the gradi- 
ent of the local magnetic flux density B perpendicular to 
the crystal surface, using the Bean model Given the 
thickness-to- width ratio of these crystals, d/w ~ 0.25, 
one has fioj s ~ 3dB/dx. 24 The crystal inhomogene- 
ity, and notably the local distribution of T c was char- 
acterize d usi ng the differential magneto-optical (DMO) 
method? 14 * 25 * Measurements in higher magnetic fields 
were performed using micron-sized Hall probe arrays,^ 
tailored in a pseudomorphic GaAlAs/GaAs heterostruc- 
ture, as well as using a Superconducting Quantum Inter- 
ference Device (SQUID)-based magnetometer. 

The vortex ensembles in several crystals, of substitu- 
tion levels x = 0.33, 0.36, and 0.49 , were imaged us- 
ing the Bitter decoration method p^ a t an applied field 
fjb H a = 20 G (2 mT) (see Section |IIIE) . In what follows, 
individual crystals will be identified as ( x = substitution 
level, sample number # ). 
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FIG. 1: (Color online) Differential magneto-optical images 
of the screening of a magnetic field /noH a = 1 G (0.1 mT) 
by BaFe2(Asi_ x P x )2 single crystals (a) x = 0.27# 2, and 
(b) x = 0.49#1. The intensity is proportional to the local 
magnetic flux density B(r). Thus, black areas are regions of 
excluded flux, while the light areas surrounding the crystal are 
traversed by the applied magnetic field. Rectangular frames 
in (a) indicate the regions where the transmittivity data of 
Fig. [5](b) are determined. The frames in (b) denote the areas 
where the decoration images of Fig. 11 (b,d) were obtained. 



III. RESULTS 
A. Spatial variation of the critical temperature T c 

Figures [I] (a,b) presents DMO images of the exclu- 
sion of an applied field fioH a = 1 G (0.1 mT) as 
one crosses the superconducting to normal transition of 
BaFe2(Asi_ iC P :E )2 single crystals x = 0.27 #2 and x = 
0.49 #1. These images reveal that T c is spatially hetero- 
geneous. While inhomogeneity is especially pronounced 
in the underdoped samples, see crystal x = 0.27 #2, it is 
also observed in the overdoped crystals. A link between 
the heterogeneity observed in underdoped samples and 
a possible phase separation between of superconducting- 
and Spin-Density Wave (SDW) antiferromagnetic phases, 
and/or that of orthorhombic and tetragonal structural 
domains is therefore not obvious. 

The T c -heterogeneity is quantified by the lo- 
cal transmittivity T H = [I(r, T) - J(r, T <C T c )\ / 
[I(r,T > T c ) - I(r,T <C T c )], extracted from the lumi- 
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FIG. 2: (Color online) (a) Tk measured on the three regions 
of crystal x — 0.27 #2 indicated in Fig.[l] (b) Transition tem- 
perature T c versus P-content x. The error bars indicate the 
local spread of T c inside a given crystal. For each x, the num- 
bering #1, #2, etc denotes different crystals from the same 
batch, (c) Dependence of the low-temperature (T = 6 K) 
critical current density on P-content x. The drawn line shows 
the evolution of the critical current due to spatial variations 
of the dopant atom density (on a scale of 5z ~ 100 nm and 
with variance Ax ~ 0.3 %) such as expected from Eq. (|9j. 



nous intensity J(r, T) in the DMO images of Fig. [T] 
T H (T) is presented in Fig. ^ for crystal x = 0.27 #2. 
The width of the superconducting transition obtained for 
the various crystals is presented in Fig. [2] b, where the 
error bar indicates the spread of T c in a given crystal, 
and the data points give the temperature where 50% of 
the crystal has become superconducting. 

B. Sustainable current density j s 

Figure [3] presents MOI of the magnetic flux pene- 
tration (after zero-field cooling) into superconducting 
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FIG. 3: Magnetic flux density distribution in 
BaFe2(Asi_ cc P cc )2 single crystals x = 0.33 #1 (a), and 
x = 0.36 #2 (b), after zero-field cooling to the indicated 
temperatures and the application of different magnetic fields. 
The top left panel of each subflgure shows the respective 
crystal, the white lines are those along which the profiles in 
Fig. [4] are extracted. The intensity in the other panels reflects 
the local flux density B(r). The bottom right-hand panel of 
(a) shows the trapped flux distribution in crystal x = 0.33 #1, 
after application and removal of /noH a = 22.1 mT. 



BaFe2(Asi_ ;E P :E )2 single crystals x = 0.33 #1, and x = 
0.36 #2, respectively. The former crystal is character- 
ized by very weak bulk pinning and, as a result, a large 
influence of geometrica l 28 * 29 * and surface barriers The 
influence of the surface screening current leads to an in- 
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FIG. 4: (Color online) Magnetic flux density profiles in 
BaFe 2 (As 1 - x P x ) 2 crystal x = 0.33 #1, at T = 19 and 26 K 
(a,b), and crystal x = 0.36 #2, at T = 19 and 22.2 K (c,d), 
after zero-field cooling and application of the magnetic field 
in successive steps AH a . The Bean-like profiles in (c,d) are 
obtained from the MOI images of Fig. [3|b). The profiles in 
(a,b) are influenced by a surface barrier and correspond to 
the crystal of Fig. [3|a). 
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FIG. 5: (Color online) Temperature dependence of the sus- 
tainable current density j s (T, B = 30 mT) , determined from 
flux density profiles obtained from MOI images, for a variety 
of crystals with substitution levels 0.23 < x < 0.49. Samples 
denoted "S" concern a limited area of the sample with the 
same number. Dashed lines are guides to the eye. 



homogeneous flux density distribution, as presented in 
Figure[3](a). In contrast, crystal x = 0.36 #2 shows regu- 
lar flux penetration, in accordance with the Bean critical 
state model. The influence of a surface barrier, present 
for both flux entry and flux exit, is also revealed by Hall- 
probe array measurements. 

The flux density profiles across the same crystals, de- 
picted in Figure [4] (a,b) and (c,d) respectively, were ex- 
tracted from the calibrated luminous intensities of the 
magneto-optical images in Fig. [3} One sees that, even 
for the same or comparable doping levels, very different 
flux density profiles can be obtained after zero-field cool- 
ing. Figure [4] (c,d) shows the Bean-like penetration of 
the magnetic flux inside crystal x — 0.36 #2, with no 
clear influence of a surface barrier, while the flux profiles 
for crystal x = 0.33 #1 in Figure [4](a,b) show, apart from 
inhomogeneity, a large discontinuity in the magnetic in- 
duction at the sample edge, characteristic of a surface 
barrier. Given the very different behavior for nearly the 
same sample composition, the origin of the bulk critical 
current density in BaFe2(Asi_ :E P :E )2 is most likely extrin- 
sic. This is supported by the temperature dependence of 
the sustainable current density j s (T, B = 30 mT) of the 
studied samples, shown in Fig. |5j The absolute value of 
j s (T) is widely disperse, even for crystals with the same 
doping level. 

In spite of the disparity, the flux pinning mechanism 
in all crystals is the same. Fig. [6] shows hysteresis loops 
of the local gradient of the magnetic induction dB/dx 
in fields of up to 2 T, obtained on crystals of different 
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FIG. 6: (Color online) Hysteresis loops of the spatial gra- 
dient dB/dx of the local magnetic induction on the surface 
of BaFe2(Asi_ cc P cc )2 single crystals of different substitution 
levels 0.27 < x < 0.58, measured using the Hall-probe mag- 
netometry technique, 26 at T — 6 K. 




\x H (T) 

FIG. 8: (Color online) (a) Magnetic hysteresis loops measured 
on BaFe 2 (Asi_ x P a ,) 2 crystal x = 0.33 #2 using a SQUID 
magnetometer, (b) Sustainable screening current j s as func- 
tion of the applied field /JoH a (in Tesla). 



composition using the Hall probe-array magnetometry 
technique.^ The hysteresis loops were measured at 6 K, 
at which flux creep has only a moderate influence. For 
all crystals, of all investigated substitution levels, one has 
the ubiquitous central peak at zero field, believed to be 
due to strong pinning by nm-scale disorder. 14 The mag- 
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FIG. 7: (Color online) Magnetic field dependence of the sus- 
tainable current density j s for crystal x = 0.36 # 2, obtained 
from the local magnetic induction gradient as measured using 
the Hall-probe magnetometry method! 2 - 6 ^ 



netic field dependence of the sustainable current density 
j s (B) was obtained from the value of dB/dx at given 
B. Fig. [7] shows j s (B) for the optimally substituted sin- 
gle crystal x = 0.36 # 2 at different temperatures. The 
j s (B) curve at the lowest T is representative of the field- 
dependent critical current density j c (B). The critical 
current density is characterized by a low-field plateau, 
followed by a j c oc B~ x l 2 decrease at higher fields, behav- 
ior that is typical o f strong flux pinning by sparse point- 
like defects. ^ 1 14 1 31 1 The contribution to j c (B) due to weak 
collective pinning of the vortex lines by atomic sized point 
pins, observed in all charge-doped iron based supercon- 
ductors as well as in Ba(Fei_ :E Ru ;E )2As2, is clearly absent 
in Fig. [7] The BaFe2(Asi_^P x )2 system can therefore be 
seen as a typical strongly pinning superconductor, in the 
sense that only large but sparse extrinsic point-like pins 
contribute to flux pinning.^ At larger temperatures and 
fields, j s (B) decreases faster than £? -1 / 2 , an effect at- 
tributed to flux creep. 

Measurements carried out to larger fields using a 
SQUID magnetometer show that, contrary to the data 
presented in Ref. [I9j the sustainable current density con- 
tinues its monotonous decrease as function of magnetic 
field - there is no "fishtail" or "second-peak" effect 14 up 
to fioH a = 5 T. The peak effect being usually associated 
with a weak collective pinning contribution to the critical 
current, we surmise that in the field range of interest the 
strong pinning mechanism is the only one at play. 
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FIG. 9: (Color online) Magnetic relaxation in 
BaFe2(Asi_ x P x )2 crystal x — 0.33 #2. (a) Relaxation 
of the magnetic flux density at three positions on the crystal 
surface - the center, and positions 20 and 40 /xm from the 
center. Data were taken at 10 K, after field cooling in 0.2 T 
and subsequently removing the field, (b) Experimental flux- 
creep activation barrier U versus j s , as obtained using the 
method outlined in Refs. 1351361 (c) Temperature dependence 
of the sustainable current density in zero applied field, and 
applied fields of 0.14, 0.2, 1, and 2 T. Measurements using 
the Hall probe-array technique (open symbols) were obtained 
from hysteresis loops such as shown in Fig. [6] while data 
from MOI are obtained from the flux-profile gradient. The 
drawn line through the j s (H a = 0)-data corresponds to the 
expected T- dependence of the depairing current; the drawn 
lines for higher fields are obtained from this by taking the 
field-dependence of the critical current ( |3|4| ) into account, 
and by correcting for flux creep using Eq. I2F. 



C. Quantitative effect of flux creep 

The influence of flux creep is assessed from relax- 
ation measurements of the local flux density using the 
Hall probe magnetometry technique. Typical exam- 
ples, shown in Fig. [9](a), show that the creep rate 
S = d\n(dB / dx) / d\nt typically amounts to a few per- 
cent. Nevertheless, j s is significantly affected by creep, 
such that it is determined as the solution of the relation 
U(j s ) = fcgTln \( tp + 1) /t], rather than by the critical 

I32I33T 



sient effects at the onset of relaxation, and r is a normal- 
ization time related to the sample inductance! 32 * 33 * The 
dependence of the flux creep barrier U (j) on current den- 
sity j can be extracted using various methods, including 
those of Maley et a/P^and Abulafia et a/.! 35 * 36 1 Applying 
the latter, we find [see Fig. |9jb)] that the creep barrier 
in optimally substituted BaFe2(Asi_ :E P :E )2 follows 



U(j) = u c 



(1) 



with values of the exponent fi ~ 1.5 — 2. Therefore, 
the time and temperature-dependence of the sustainable 
screening current density is described by 



Jc 



In 



t + t 



(2) 



The impact of flux creep on the temperature dependence 
of the sustainable current is depicted in Fig.[9jc), which 
shows j s (T)-curves for crystal x = 0.33 #2, for dif- 
ferent B. The curve in zero applied field is little af- 
fected by creep, and roughly follows the expected tem- 
perature dependence of the depairing current, j(0,T) ~ 
e (T)/$o£(T) ($0 = h/2e is the flux quantum). Here, 
the vortex line energy, So(T) = <£o/47r/ioA 2 6 , propor- 
tional to the superfluid density n s ~ A~ 2 , is evalu- 
ated using the data for the in-plane penetration depth 
X ab (T) of Ref. SU and the coherence length f(T) ~ 
^(0)v / ( 1 + T / T c)/(1 -T/T c ). The curves for varying 
applied field can then be well described by taking the 
creep barrier prefactor U C (T) oc €o(T), j c (B,T) oc 
j(0, T)£? -1 / 2 , fj, = 1.6, and In [(t + t ) /r] = 20P There- 
fore, the temperature dependence of the screening cur- 
rent in fields larger than 0.1 T is essentially determined 
by flux creep. 



D. Extraction of pinning parameters 

We now analyze the j s (B)-curves measured at low 
temperature, representative of the critical current density 
j c (B), and which bear the hallmarks of str ong p inning. 
These are the plateau at low magnetic 

followed by a power-law decrease as a function of the flux 
density ffj 14 l 31 l which can be described as 



MB) 



fp £\ao 



current density j j 1 'Here, to is a time describing tran- 



1/2 

The main parameter, C = {irU p ni/ S\Sq) 1 , is the av- 
erage distance between effective defects pinning a single 
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FIG. 10: (Color online) Average distance C between effective 
pins versus P content. C is obtained from j s (B)-data at 5 K 
(red bullets •), and from Bitter decoration at Tf = 0.87T C 
(blue squares ■). The drawn lines show the comparison with 
Eq. (10), for variations of the dopant atom density Ax = 
0.3 %, on a characteristic length scale Sz = 100 nm. 
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FIG. 11: Bitter decoration images of BaFe2(Asi_ cc P cc )2 single 
crystals for an applied field /noH a = 20 G (a,c) x = 0.36 #2, 
and (b,d) x = 0.49 #1. The white arrows indicate the vortex- 
free Meissner belt observed (a) near a surface step and (d) 
at the edge of the crystal x = 0.49 #1. The regions where 
images (b) and (d) were obtained with respect to the entire 
crystal are depicted in Fig. [T] 



E. Vortex imaging by Bitter decoration 



vortex in the low-field limit The crossover field B* 
is that above which the intervortex repulsion limits the 
number of effective pins per vortex, and f p is the maxi- 
mum pinning force exerted by a single strong pin. U P /[J] 
is the pinning energy of a single strong pin, rii is the pin 
density, ao = (Qo/B) 1 / 2 is the intervortex distance, and 
£\ = \ a b j X c is the penetration depth anisotropy. 



Equations 
ning force 

-3/2 



(3) 
of 



and Q 
a single 



show that 
strong pin 



the 

fp 



pin- 



(% /2 exM{j 2 (0)/[dj(B)/dB- 1 / 2 ]} can be obtained 
from the experimentally measured low temperature, low- 
field current density j c (0), and the slope dj c (B)/dB~ 1 / 2 
at intermediate fields. Reserving our attention to the 
crystals used in the Bitter decoration experiments pre- 
sented below, we obtain, for an estimated £\ = 0.15p^ 
f p ^8x 10~ 13 N for both crystals (x = 0.36#2) and 
(x = 0.49#1). This value is twice larger than that 
measured in Ba(Fei_ £C 00^)2 AS2P An evaluation at 
the highest measurement temperature of 0.8T C yields 
f p = 2 x 10 -14 N; however, this value is likely to be 
overestimated due to creep. Similarly, one can extract 

shows that C is 
in ac- 



10 



the length C = f p /$oj c {0). Figure 
of the order of several dozen to hundreds of nm 
cordance with the strong pinning hypothesis. Moreover, 
the distance between effective pins clearly increases as 
function of P-content x. 



The Bitter decoration techniqu e! 2 1 27 1 was used to image 
the vortex ensemble in BaFe2(Asi_ cc P cc )2 single crystals 
with three different substitution levels, x = 0.33, x = 
0.36, and x = 0.49. The experiments were realized under 
field cooled (FC) conditions, with a field fioH a = 20 G 
(2 mT) applied parallel to the c-axis of the crystals. The 
decoration experiment for the crystal with x — 0.33 was 
not successful, presumably due to the large value of the 
penetration depth.^ 

The vortex configurations shown in Fig. [TT] for crystals 
x = 0.36 #2, and x = 0.49 #1, are representative of what 
is observed over the entire specimens. From the decora- 
tion images, we obtain, for both crystals, the average 
value of the magnetic induction as B int = n v $>o w 19 G 
(with n v the vortex density). This is 1 G smaller than 
the applied field during the experiment. Notwithstand- 
ing Ref. [39j there is therefore evidence for Meissner ex- 
clusion of the magnetic flux. Moreover, the Meissner cur- 
rent manifests itself as a vortex- free "Meissner belt" along 
the edges of decorated crystal x = 0.49 #1, as well as 
near surface steps that appear during preliminary cleav- 
age of the samples, as indicated in Figure [TT] Long vortex 
chains reminiscent of those observed in Ref. [TT] are also 
observed in the decoration images, for both investigated 
P-contents, however the chains are more pronounced in 
crystal x = 0.36 #2. 

Figures [l2] (a) and (b) present the respective Delaunay 
triangulations of the vortex ensembles of Fig. 11 (a) and 
(b). Here, blue dots represent vortices with sixfold coor- 
dination, while red dots represent vortices with different 
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FIG. 12: (Color online) Delaunay triangulation of vortex en- 
sembles observed in BaFe2(Asi_ cc P cc )2 single crystals (a) x = 
0.36 #2 and (b) x = 0.49 #1, and presented in Figs. 11 (a,c) 
respectively. Blue dots represent vortices with sixfold coordi- 
nation, while red dots represent differently coordinated vor- 
tices. The insets show the Fourier transform of the vortex 
positions. (c,d) Nearest neighbor distance distributions for 
the respective triangulations. The insets present the respec- 
tive coordination number distributions. 
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FIG. 13: (Color online) (a-d) Color-coded maps of the nor- 
malized individual v orte x interaction energies calculated from 
the images of Fig. ll'a-d) using Eq. |5]), and represented 
in the same configuration. Lower panels: normalized in- 
teraction energy distributions for BaFe2(Asi_ cc P cc )2 crystals 
(x = 0.36 #2) (e,g), and (x = 0.49 #1) (f,h). Vortices located 
within a distance 10A a b from the map edge are excluded from 
the histograms. The interaction energy per vortex of the tri- 
angular lattice (^-function) is represented by the central or- 
ange line in each histogram. 



coordination number. Even if there are some regions with 
several adjacent sixfold coordinated vortices, the insets to 
Fig. 12 (c) and (d) reveal that the latter represent less 
than half of the total (46% for crystal x = 0.49 #1, and 
43% for crystal x = 0.36 #2). In fact, for the magnetic 
field of 20 G under study, the percentage of sixfold coor- 
dinated vortices is the same as in Ba(Feo925Coo. 075)2^82 
(43%) 40 and in Ba .6Ko.4Fe 2 As2 (44%). 11 Moreover, the 
coordination number histograms have the same width for 
the three materials, revealing similar disorder of the vor- 
tex ensemble. The absence of vortex lattice order is fur- 
ther brought out by the Fourier transforms of the vor- 



tex positions shown in the insert to Fig. [l2] (a) and (b). 
Nevertheless, the spatial distribution of vortices in the 
two BaFe2(Asi_ :c P :c )2 crystals presents smaller density 
fluctuations than previously observed in the Co-doped 
materialP Panels (c,d) of Fig. 
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show the distributions 
of nearest-neighbor intervortex distances. These have 
a mean value |r^ | = 1 /im, while the lattice parame- 
ter for a tria ngular perfect lattice of the same density is 
a A = 1.075 y/$ /B = 1.12 /im. This shift is due to the 
existence of more densely packed vortices, notably in the 
chain-like structures. 
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FIG. 14: (Color online) Normalized color-coded maps of the 
modulus of the pinning force (per unit length), calculated 
from the images of Fig. 11 using Eq. Lower panels: pin- 
ning force distributions for BaFe2(Asi_ a; P :c )2 single crystals 
[x = 0.36 #2), (e,g), and (x = 0.49 #1), (f,h). Again, vor- 
tices located within a distance 10X a b from the map edge are 
excluded from the histograms. 



IV. DISCUSSION 

In what follows, we adopt the procedure of Ref.[2]to de- 
termine the vortex interaction energy. For this, one needs 
to know the value of A a 5 at the temperature Tf at which 
the vortex ensemble is frozen. As in Ref. 2, we use the in- 
formation that can be obtained from vortex lines situated 
near surface steps. Such steps may act as barriers, but, 
due to the circulation of the Meissner current, they also 
prevent vortex lines from being situated right at their 
edge. Inserting the height of the surface step in Fig. [TT] 
(a), h = 1.3 /im, and the width of the vortex-free region 
close to the step, u = 1.2 /im, in Eq. (1) of Ref. O the 
value of the penetration depth at the freezing tempera- 
ture is graphically estimated as A a 5(Tj) w 700 nm. Using 
the temperature dependence of A a b(T) from Ref. |4T| one 



obtains the freezing temperature of the vortex ensemble 
as Tf « 0.87T C . Thus, even though vortices are frozen at 
a relatively high reduced temperature, its value is lower 
than T f /T c = 0.95 found in Ba(Fei_ x Co :c )2As2. 

Using \ a b(Tf) and the vortex positions extracted from 
Fig. [TT] the interaction energies of the individual vortices 
can be calculated as 



(5) 



Here Kq(x) is the lowest-order modified Bessel function, 
and \rij\ is the distance from vortex i to vortex j. For 
each vortex z, only neighbors j contained within a circle 
of radius of 10A a &(T) are taken into account. This radius 
was chosen after verification that vortices situated at a 
larger distance do not significantly contribute to £\ nt . 

Figures [T3^a-d) present the vortex interaction energies 
as color-coded maps, with the energy scale normalized 
by eo(Tf), as extracted from the images of Fig. 
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to (d) respectively. The maps show a globally homo- 
geneous distribution; however, a number of denser re- 
gions exist. Histograms of the interaction energies for 
the maps (a) to (d) are presented in Fig. 13 (e) to (h). 
Note that the presence of the chain-like features with 
a denser vortex arrangement broadens the histograms 
for BaFe 2 (As 1 _ x P x ) 2 single crystal (x = 0.36 #2). Still, 
the energy distributions are considerably narrower than 
those found in Ba(Fei -^Co^) 2 AS2P Furthermore, all dis- 
tributions are centered about the average Sint *~ 3.5so? 
which corresponds to the interaction energy value (5- 
peak) of the perfectly traingular Abrikosov lattice for 
this particular vortex density. Therefore, in contrast to 
Ba(Fei_ ;E Co ;E )2As2P no pinning-induced shift of the av- 
erage value of the energy distribution histogram with re- 
spect to the 5-peak value is observed. 

The large shift found in Ref. 2 was interpreted in 
terms of a large average pinning energy in the vicinity 
of T c , which can only be the result of T c heterogene- 
ity. The absence of such a shift in isovalently substi- 
tuted BaFe2(Asi_ ;E P ;E )2 suggests that spatial inhomo- 
geneity of T c is irrelevant for vortex pinning in this ma- 
terial. A probable reason for this is a smoother temper- 
ature dependence of the pinning potential. The pinning 
energy due to spatial variations of the vortex line en- 
ergy, on a scale Sz and with variance Asq, writes U p ~ 
Ae (T)Sz - [An s (0)/n s (0)][de (T)/dn s (0)]n s (T) + 
[e /n s (0)}[dn s (T)/dT c }AT c . At temperatures close to T c , 
the second contribution, due to spatial fluctuations of T c , 
dominates the pinning energyP However, whereas in ma- 
terials such as Ba(Fei_ a; Co :c )2As2 - in which n s is nearly 
linear in 1 — T/T c as one approaches T c (see Fig. 15)- 



this contribution is large, in BaFe2(Asi_ :E P :c )2 with a 
smoother temperature dependence of n s this contribu- 
tion vanishes. In fact, the T c heterogeneity within the 
decorated areas of the P-substituted material, such as 
observed by the DMO technique, does not result in qual- 
itatively different vortex arrangements in different parts 
of the crystal [see Figure [T](b)]. 
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FIG. 15: Temperature dependence of the superfluid den- 
sity, n s (T)/n s (0) oc A^ b (0)/A^(T) for various doping levels 
x of Ba(Fei_ x Co x )2As2 (upper curves, small data points), of 
Ba .55Ko.45Fe 2 As2 (•) and BaFe 2 (As .67Po. 33)2 ( ) (both from 
Ref.HU). 



Figure [14] shows maps of the modulus of the pinning 
force for each individual vortex, calculated from 



2£Q Tjj 
^ab \ T ij\ 



following a similar procedure as used for the determi- 
nation of E\ nt . Here K\(x) is the first order modified 
Bessel function. Since the rendered vortex configurations 
in Fig. [TT] are in a stationary state at the freezing tem- 
perature Tf, the calculated intervortex repulsive force 
must be balanced by the pinning force. The maps of 
Fig. [l4^a) to (d) therefore represent the minimum local 
pinning force for each vortex, min(|f^|). The distribu- 
tions of min(|fj|) shown in Fig. [l4|(e) to (h) allow one 
to estimate the average pinning force per vortex and per 
unit length. We obtain |f^| ~ 3.5 x 10 -6 Nm" 1 for crys- 
tal x = 0.49 #1 and |f;| - 4.5 x 1(T 6 Nm" 1 for crystal 
x = 0.36 #2. These (high temperature) values are com- 
parable to those found in Ba(Fei_a; 00^)2 AS2P 

The maps of the local pinning force moduli and the 



interaction energy shown in Figs. [13] and [14| respectively 
are correlated. However, the respective probability dis- 
tributions are clearly broader for the lower substitution 
level x = 0.36. As in Ref. [2j the ratio of the elementary 
pinning f orce p er pin f p: extracted from the j(B) curves 
and the value of I 



in section 



III A 



\fi \ obtained from Bitter 
decoration, allows one to evaluate an upper bound on C 
in an independent manner. Using the low-temperature 
value f p ^Sx 10 -13 N yields C = 180 nm for crystal 
(x = 0.36 #2) and I = 230 nm for crystal (x = 0.49 #1). 
Figure [10] shows that these numbers are consistent with 
those directly extracted from the sustainable current den- 
sity measurements. 

Several reasons can be invoked to explain the observed 



enhancement of C with increasing P-content x. First is 
the increase of the superfluid density n s (x) as function 
of xW^ Other possibilities are a decrease of the density 
rti of pinning centers for larger x, and the decrease of the 
penetration depth anisotropy e^ 1 for higher substitution 
levels. We first investigate the effect of the increase of 
the superfluid density with x. Assuming that nm-scale 
fluctuations of the dopant atom density (with variance 
Ax) are spatially isotropic, their effect on the pinning 
force 



</ P >~</ Ve (r)dz) ~ Ae 

J5z 



(7) 



and the pinning energy U p ~ AsqSz can be estimated by 
exploiting the dependence \ a b(%)W^ 



(f P ) ~ As 



de dX ab 



Ax. 



dXab dx 

The low-field critical current density Q becomes 



(8) 



1/2 



de d\ ah 



dXab dx 



Ax 



3/2 



(9) 



(6) the length between effective pinning centers 



£0 



7T 1 / 2 \riiAeoSz 



2tt 1 /2 (mSz) 1 / 2 



1/2 



1 de dXat 



e d\ ab dx 



Ax 



-1/2 



. (10) 



Assuming that for density fluctuations of the dopant 
atoms the pin density rii ~ (5z)~ 3 scales as the inverse 
cube of the defect size, these expressions can be directly 
compared to the dependence of the critical current den- 
sity on P-content (see Fig. [2]:), as well as that of C (see 
Fig. [TO]). For Sz = 100 nm and Ax = 0.3 %, the qualita- 



tive trend with x of both quantities can be reproduced. 
Thus, the dependence of the penetration depth on P- 
content accounts, at least in the overdoped regime, for an 
enhanced probability of encountering larger critical cur- 
rent densities around optimal substitution. The similar 
dispersion of critical current density and pinning length 
data for different x around the model estimations ([9| and 
( [To] ) in Figs. [2 2 and 10 would then be due to similar disor- 
der of the P-distribution in the different growth batches, 
eliminating the need to invoke a decreasing disorder with 
increasing x. Still, the evolution C(x) is reminiscent of 
that of the normal-state mean free-path, extracted by 
Shishido et. at. from de Haas-van Alphen oscillations 
of the Landau magnetization!^ They reported that the 
mean free-path for the (3 orbits in BaFe2(Asi_ x P iC )2 sin- 
gle crystals increases from I ~ 170 A to 800 A when P 
content varies from x = 0.41 to 1. This would imply the 
presence of structural defects that act as strong pinning 
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centers that are unrelated to spatial composition fluctu- 
ations. Last but not least, the increase of e\(x), as in- 
ferred from the doping-dependent evolution of the Fermi 
surface,^ and corresponding to a smaller Fe-pnictogen 
distance for large x, would lead to stiff er vortex lines and 
to less pinning for higher P-content. 

Finally, we address the clear presence of chain-like 
structures in the Bitter decoration images. The following 
hypotheses can be suggested to explain the appearance 
of these chains. The first is the possibility of native het- 
erogeneity of the crystals under study, introduced during 
growth. The presence of line-like defects, or of linear 
agglomerates of point-like defects giving stronger local 
pinning, will lead to the appearance of vortex alignments 
or chainlike structures, much as this was found in e.g. 
Bi2Sr2CaCu208+<$P^"^ A second possibility is that, the 
vortex images having been obtained after field-cooling, 
and since the BaFe2(Asi_ :c P :E )2 material under consider- 
ation show Meissner expulsion, a certain fraction of vor- 
tices must exit the material before the vortex ensemble 
is frozen at Tf. The chains possibly correspond to flow 
channels for these exiting vortices, the flux in interme- 
diate areas remaining pinned!^ Finally, the multi-band 
character of BaFe2(Asi_ iC P :c )2 may be responsible for the 
occurrence of vortex chains, with the possible existence 
of an attractive part in the inter- vortex potential!^ We 
leave these questions open for further work. 

V. SUMMARY AND CONCLUSION 

We have presented an overview of vortex pinning in 
single crystals of the isovalently substituted iron-based 
superconductor BaFe2(Asi_ :c P :E )2, in which we have at- 
tempted to correlate the sustainable screening current 
density as a function of temperature, field, and doping 
x, with the structural properties of the vortex ensemble. 
The critical current density in BaFe2(Asi_ ;E P ;E )2 is, over- 
all, very well described by the strong pinning scenario of 
Ref. [3TJ which allows one to extract elementary pinning 
forces (of the order of 10 -13 N) and the distance between 
effective pins. The latter is of the order of 100 nm, and 
increases as a function of doping level x. These values are 
consistent with those independently obtained by means of 
the magnetic decoration technique. Contrary to Ref. [19j 
we find no contribution of weak collective pinning to the 
sustainable current density, suggesting that P-atoms are 
not responsible for quasi-particle scattering.^ The sus- 



tainable current data are affected by flux creep, which 
prohibits one from drawing definite conclusions concern- 
ing the temperature dependence of pinning. 

Bitter decoration reveals slightly more ordered vor- 
tex ensembles (the number of sixfold coordinated vor- 
tices is slightly higher) than those observed in charge- 
doped Ba(Fei_ iC Co :E )2As2. Also, the interaction energy 
and pinning force distributions in BaFe2(Asi_ x P :c )2 are 
much narrower than those in Ba(Fei_ :E Co ;E )2As2, and are 
not shifted with respect to the interaction energy of a 
perfectly triangular vortex lattice with the same density. 
These observations exclude a role of spatial variations 
of the critical temperature T c in determining the frozen 
vortex state obtained upon field-cooling. The absence 
of the weak- collective pinning contribution to the criti- 
cal current density in BaFe2(Asi_ x P :E )2 means that the 
strong pinning contribution is what generates the disor- 
dered vortex configurations. 

The main features of strong vortex pinning in 
BaFe2(Asi_ cc P :c )2, such as the energy- and force his- 
tograms, the density of effective pins, and the evolution 
of the critical current density with P-content, can be un- 
derstood using a model for heterogeneity of the superfluid 
density on the scale of several dozen to several hundred 
nanometers, due to an inhomogeneous distribution of the 
dopant atoms. A small spatial variance Ax ( of the or- 
der of 0.3 % ) suffices to explain the magnitude of j c . 
This model explains why pinning is of extrinsic origin, 
and why the disorder can be readily annealed. 12 Even 
if the x-dependence of the critical current density can 
be well described without any assumption of less disor- 
der in crystals with higher P-content, the reminiscence 
of the evolution with P-content of the mean distance be- 
tween effective pinning sites and the normal-state mean- 
free patrP^ suggests that a second type of pinning centers 
may be at play. This could be areas of enhanced local 
strain, more prominent at optimal doping, such as these 
arise from the very different Fe-As and Fe-P distances.^ 
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